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Abstract 



We construct and discuss a "realistic" example of SUSY SU{5) GUT model, with an addi- 
tional U{1) flavour symmetry, that is not plagued by the need of large fine tunings, like those 
associated with doublet-triplet splitting in the minimal model, and that leads to an accept- 
able phenomenology. This includes coupling unification with a value of as{mz) in much better 
agreement with the data than in the minimal version, an acceptable hierarchical pattern for 
fermion masses and mixing angles, also including neutrino masses and mixings, and a proton 
decay rate compatible with present limits (but the discovery of proton decay should be within 
reach of the next generation of experiments). In the neutrino sector the preferred solution is 
one with nearly maximal mixing both for atmospheric and solar neutrinos. 
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1 Introduction 



The idea that all particle interactions merge into a unified theory at very high energies is 
so attractive that this concept has become widely accepted by now. The quantitative success 
of coupling unification in Supersymmetric (SUSY) Grand Unified Theories (GUT's) has added 
much support to this idea |^. The recent developments on neutrino oscillations pointing 
to lepton number violation at large scales, have further strengthened the general confidence. 
However the actual realization of this idea is not precisely defined. On the one hand, at the 
Planck scale, Mp^, the unification of gravity with gauge interactions is a general property of 
superstring theories. On the other hand, in simple GUT models gauge interactions are unified 
at a distinctly lower mass scale Mqut within the context of a renormalizable gauge theory. But 
this neat separation between gauge unification and merging with gravity is not at all granted. 
The gap between Mqut and Mpi could be filled up by a number of threshold effects and 
several layers of additional states. Also, coupling unification could be realized without gauge 
unification, as suggested in some versions of superstring theory or in fiipped SU{b). 

Assuming gauge unification, minimal models of GUT's based on S'[/(5), 5*0(10), ... have 
been considered in detail. Also, many articles have addressed particular aspects of GUT's 
models like proton decay, fermion masses and, recently, neutrino masses and mixings. But 
minimal models are not plausible as they need a large amount of fine tuning and are therefore 
highly unnatural (for example with respect to the doublet-triplet splitting problem or proton 
decay). Also, analyses of particular aspects of GUT's often leave aside the problem of embed- 
ding the sector under discussion into a consistent whole. So the problem arises of going beyond 
minimal toy models by formulating sufficiently realistic, not unnecessarily complicated, rela- 
tively complete models that can serve as benchmarks to be compared with experiment. More 
appropriately, instead of "realistic" we should say "not grossly unrealistic" because it is clear 
that many important details cannot be sufficiently controlled and assumptions must be made. 
The model we aim at should not rely on large fine tunings and must lead to an acceptable 
phenomenology. This includes coupling unification with an acceptable value of a;s(m^), given 
a and sin^Ow at m^, compatibility with the more and more stringent bounds on proton decay 

1^, agreement with the observed fermion mass spectrum, also considering neutrino masses 
and mixings and so on. The success or failure of the program of constructing realistic models 
can decide whether or not a stage of gauge unification is a likely possibility. 

Prompted by recent neutrino oscillation data some new studies on realistic GUT models 
have appeared in the context of 5*0(10) or larger groups 0. In the present paper we address the 
question whether the smallest SUSY SU (5) symmetry group can still be considered as a basis 
for a realistic GUT model. We indeed present an explicit example of a realistic SU (5) model, 
which uses a U{1) flavour symmetry as a crucial ingredient. In principle the flavour symmetry 
could be either global or local. We tentatively assume here that the flavour symmetry is global. 
This is more in the spirit of GUT's in the sense that all gauge symmetries are unifled. The 
associated Goldstone boson receives a mass from the anomaly. Such a pseudo-Goldstone boson 
can be phenomenologically acceptable in view of the existing limits on axion-like particles 0. 

In this model the doublet-triplet splitting problem is solved by the missing partner mech- 
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anism stabilized by the flavour symmetry against the occurrence of doublet mass lifting 
due to non renormalizable operators. Relatively large representations (50, 50, 75) have to be 
introduced for this purpose. A good effect of this proliferation of states is that the value of 
0-s{f^z) obtained from coupling unification in the next to the leading order perturbative ap- 
proximation receives important negative corrections from threshold effects near the GUT scale. 
As a result, the central value changes from ^^(m^) 0.130 in minimal SUSY SU{b) down to 
Ois{fTiz) ~ 0.116, in better agreement with observation 0, |Tl[. The same U{1) flavour sym- 
metry that stabilizes the missing partner mechanism is used to explain the hierarchical structure 
of fermion masses. In the neutrino sector, the mass matrices already proposed in a previous 
paper by two of us are reproduced [|12| . The large atmospheric neutrino mixing is due to a large 



left-handed mixing in the lepton sector that corresponds to a large right-handed mixing in the 
down quark sector. In the present particular version maximal mixing also for solar neutrinos 
is preferred. A possibly problematic feature of the model is that, beyond the unification point, 
when all the states participate in the running, the asymptotic freedom of SU{b) is destroyed 
by the large number of matter fields. As a consequence, the coupling increases very fast and 
the theory becomes non perturbative below Mpi. In the past models similar to ours have been 
considered, but were discarded just because they contain many additional states and tend to 
become non perturbative between Mqut and Mpi. We instead argue that these features are not 
necessarily bad. While the predictivity of the theory is reduced because of non renormalizable 
operators that are only suppressed by powers of Mqut/ ^ with A < Mp;, still these corrections 
could explain the distortions of the mass spectrum with respect to the minimal model, the 
suppression of proton decay and so on. However, it is certainly true that also in this case, as for 
any other known realistic model, the resulting construction is considerably more complicated 
than in the corresponding minimal model. 



2 The Model 

The symmetry of the model is SUSY SU (5) ®U{1). We call Q the charge associated with U{\). 
The superpotential of the model has three parts: 

W = W1+W2 + WZ . (1) 

The Wi term only contains the field Y of the representation 75 of SU{^) with Q = Q: 

wi = ci + My Y'^ . (2) 

The effect of w\ is to provide Y with a vev of order My /ci My ~ Mqut and to give a mass 
to all physical components of Y , i.e. those that are not absorbed by the Higgs mechanism (the 
75 uniquely breaks SU{h) down to SU{?>) ® SU{2) (g) U{1)). 

The W2 term induces the doublet-triplet splitting: 

W2 = C2 HYH,o + csHYH^ + C4 H,oH^X . (3) 
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Here H and H are the usual pentaplets of Higgs fields, except that now they carry non opposite 
Q charges. It is not restrictive to take C2, C3 and C4 real and positive. The renormalizable 
couplings that appear in W2 are the most general allowed by the SU{5) and Q assignments of 
the fields in eq. which are given as follows: 



field 


Y H 


H 


H50 


H50 


X 


SU{5) 


75 5 


5 


50 


50 


1 


Q 


-q 


q-l 


Q 


1-q 


- 1 



(4) 

The value of q will be specified later. At the minimum of the potential, in the limit of un- 
broken SUSY, the vevs of the fields H, H, H^q and H-^ all vanish, while the X vev remains 
undetermined. As we shall see, when SUSY is softly broken the light doublets in H and H 
acquire a small vev while the X vev will be fixed near the cut-off A, of the order of the scale 
between Mqut and Mpi where the theory becomes strongly interacting (we shall see that we 
estimate this scale at around 10 — 20 Mqut, large enough that the approximation of neglect- 
ing terms of order Mgut/^ is not unreasonable). The missing partner mechanism to solve 
the doublet-triplet splitting problem occurs because the 50 contains a (3, 1), i.e. a coloured an- 
titriplet, SU{2) singlet (of electric charge 1/3) but no colourless doublet (1,2). The U{1) flavour 
symmetry protects the doublet Higgs to take mass from radiative corrections because no HH 
mass term is allowed. Also no non renormalizable terms of the form HHY"^X"^ (m, n > 0) 
are possible, because X has a negative Q charge. This version of the missing partner mech- 
anism was introduced in ref. and overcomes the observation in ref. that, in general, 
non renormalizable interactions spoil the mechanism. The Higgs colour triplets mix with the 
analogous states in the 50 and the resulting mass matrix is of the see-saw form: 



TTLt 



C2{Y) 
C3{Y) c,{X) 



Defining = C4(X) the eigenvalues of the matrix mTTfir are the squares of: 

1 



mTi,2 = - 



ml + (C2 + C3)2(F)2 ± Jml + {c^ - cs^iY) 



(5) 



(6) 



Note that mTi'mT2 = C2C3{Y)'^. The effective mass that enters in the dimension 5 operators 
with \AB = 1| is 

mTimT2 C2C3 (F)^ 



rriT 



C4 {X) 



(7) 



50 



The W3 term contains the Yukawa interactions of the quark and lepton fields \l'io, and 
^'i, transforming as the representations 10, 5 and 1 of SU(5) respectively. We assume an exact 
i?-parity discrete symmetry under which \&io, ^E^g and \E'i are odd whereas H, H, H^q and H-, 
are even. The term is symbolically given by 

W3 = ^ioGu{X,Y)^ioH + ^ioGd(X,F)^5i? + '^-,G,{X,Y)^,H 
+ M^iG'A/(X,r)^i + ^,oG^iX,Y)^^oH^ . 
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The Yukawa matrices G^, Gd, G^, Gm and G-^ depend on X and Y and the associated mass 
matrices on their vevs. The last term does not contribute to the mass matrices because of 
the vanishing vev of H-^, but is important for proton decay. The pattern of fermion masses is 
determined by the U{1) flavour symmetry that fixes the powers of A = {X)/A for each entry 
of the mass matrices. In fact X is the only field with non vanishing Q that takes a vev. The 
powers of A in the mass terms are fixed by the Q charges of the matter field \E' and of the Higgs 
fields H and H. We can then specify the charge q that appears in and the Q charges of the 
matter fields \E' in order to obtain realistic textures for the fermion masses. We choose q = 2, 
so that we have from the table in fHI): 



Q{H) = -2 and Q{H) = 1 

and, for matter fields 

g(^io) = (4,3,l) , Q(vi/5) = (4,2,2) , 
The Yukawa mass matrices are of the form: 

Gr{{X), = X-^^Grmk, 



:i,-i,o) 



M, d, u, M 



(9) 



(10) 



(11) 



We expand Gr{{Y))ij in powers of (Y) and consider the lowest order term at first. Taking 
Gr{0)ij of order 1 and riij as dictated by the above charge assignments we obtain [|: 



rriu 



1 

V2 



1 
72 



"A^ 


A^ 


A3 


A^ 


A^ 


A2 


_A3 


A2 


1 


"A3 


A 


A2" 


A 





1 


A 





1 



rrid 



1 

V2 



m 



ma] 





A3 


A3" 






X' 


A^ 


A2 


VdX"^ , 


x^ 


1 


1 _ 








"A2 


1 


A" 






1 








M . 




A 





1_ 




need A ~ 


Xc 




0.22, Xc 



(12) 
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being 

with two important 
The factor 



g the 



the Cabibbo angle. These mass matrices closely match those of ref. 
special features. First, we have here that tan/3 = Vu/vd ~ rrit/mbX'^, which is small. 
A"^ is obtained as a consequence of the Higgs and matter fields charges Q, while in ref. 
H and H charges were taken as zero. We recall that a value of tan/3 near 1 is an advantage 
for suppressing proton decay. A small range of tan j3 around one is currently disfavored by the 
negative results of the SUSY Higgs search at LEP [|15|]. Of course we could easily avoid this 
range, if necessary. Second, the zero entries in the mass matrices of the neutrino sector occur 
because the negatively Q-charged X field has no counterpart with positive Q-charge. Neglected 
small effects could partially fill up the zeroes. As explained in ref. |]T^ these zeroes lead to 
near maximal mixing also for solar neutrinos. 



md 



A problematic aspect of this zeroth order approximation to the mass matrices is the relation 
This equality is good as an order of magnitude relation because relates large left- 



m 



T 



handed mixings for leptons to large right handed mixings for down quarks [|T^, |T6[|. However 
^ In our convention Dirac mass terms are given by Li^raR* and the light neutrinos effective mass matrix is 



4 



the implied equalities mf,/mT- = rris/m^ = md/me = 1 are good only for the third generation 
while need to be corrected by factors of 3 for the first two generations. The necessary corrective 
terms can arise from the neglected terms in the expansion in (Y) of Gr{{Y))ij |T^. The higher 



order terms correspond to non renormalizable operators with the insertion of n factors of the 
75, which break the transposition relation between rrid and rrie. For this purpose we would like 
the expansion parameter {Y)/A to be not too small in order to naturally provide the required 
factors of 3. We will present in the following explicit examples of parameter choices that lead 
to a realistic spectrum without unacceptable fine tuning. 

The breaking of SUSY fixes a large vev for the field X by removing the corresponding flat 
direction, gives masses to s-partners, provides a small mass to the Higgs doublet and introduces 
a n term. Up to coefficients of order one we can write down the terms that break SUSY softly: 



soft 



2| |2 I 2| |2 

m \x\ + m c 



+ 



2|-|2 

m c 



+ m(ccx + h.c.) + 



(13) 



where c, c and x denote the scalar components of H^q, H-^ and X respectively, and dots stand 
for the remaining soft breaking terms, including mass terms for the scalar components of the 
Higgs doublet fields 0. In the SUSY limit and neglecting the mixing between the (50, 50) and 
the (5,5) sectors, fermions and scalars in the 50 and 50 have a common squared mass c\x'^ . 
When SUSY is broken by the soft terms for each fermion of mass c^x there are two bosons of 
squared masses c\x'^ + m? ± mx. The x terms in the scalar potential at one loop accuracy are 
given by: 



V = m^x^ + 



50 
64^ 



2(c'ix^ + + mx)^ 



log 



clx"^ + + mx ' 



+ 2(caX +m — mx) 



log 



-^2 ^2 




clx"^ + — mx" 
A2 



A2 
3/2 



3/2 



(14) 



A numerical study of this potential in the limit m <^ A and for C4 of order one, shows that 
the minimum is at x of order A (somewhat smaller than A but close to it). The expression 
for V in eq. (|T^) provides a good approximation of the scalar potential only in a small region 
around A. Outside this region the perturbative approximation will break down. We take our 
numerical analysis as an indication that the minimum occurs near A and we assume that the 
true minimum occurs at x = {X) = 0.25A, that is A = 0.25. The complex field x describes two 
physical scalar particles. The one associated to the real part of x has a mass of order m and 
couplings to the ordinary fermions suppressed by 1/A. The particle associated to the imaginary 
part of X is massless only in the tree approximation. Due to the anomaly of the related U{1)q 
current the particle acquires a mass of order f.„mT,/{X). Cosmological bounds on (X) have 
been recently reconsidered in ref. where it has been observed that {X) of the order of the 
grand unification scale is not in confiict with observational data. 



^To find the minima of the scalar potential it is not restrictive to set to zero the imaginary part of x, which 
will be understood in the remaining part of this section. 



5 



An alternative possibility is to assume that the U{1) symmetry is local |jT3[. In this case 
the super symmetric action contains a Fayet-Iliopoulos term and the associated D-term in the 
scalar potential provides a large vev for x, of the order of the cut-off scale A. 

A fi term for the fields H and H of the appropriate order of magnitude can be generated 
according to the Giudice-Masiero mechanism []TB[. Assume that the breaking of SUSY is induced 
by the 9'^ component of a chiral (effective) superfield S, singlet under SU{5) with Q = 0. In 
the Kahler potential a term of the form 

K = +h.c. (15) 

is allowed. The vevs of S and X, {S) ~ 6'^mMpi, {X) = AA lead to an equivalent term in the 
superpotential of the desired form fiHH with ji ~ XmMpi/A which can be considered of the 
right order. 



3 Coupling Unification 

It is well known that in the minimal version of SUSY SU(5) the central value of as{mz) required 
by the constraint of coupling unification is somewhat large: as{mz) ~ 0.13 |P, |10|, |TT|, |19|, pO| . 
In the model discussed here, where the doublet triplet splitting problem is solved by introducing 
the SU{5) representations 50, 50 and 75, the central value of as{mz) is modified by threshold 
corrections near Mqut which bring the central value down by a substantial amount so that 
the final result can be in much better agreement with the observed value. As discussed in refs. 

0, this remarkable result arises because the 24 of the minimal model is replaced by the 
75. The mass splittings inside these representations are dictated by the group embedding of 
the singlet under SU{?>) ® SU{2) ® U{1) that breaks SU(5). The difference in the threshold 
contributions from the 24 and the 75 has the right sign and amount to bring ^^(m^) down even 
below the observed value. The right value can then be obtained by moving ttisusy and uit 
in a reasonable range. The difference in the favored value of rriT with respect to the minimal 
model in order to reproduce the observed value of as{mz) goes in the right direction to also 
considerably alleviate the potential problems from the bounds on proton decay. Note that there 
are no additional threshold corrections from the 50 and 50 representations because the mass of 
these states arise from the coupling to the field X which is an SU{5) singlet. Thus there are 
no mass splittings inside the 50 and the 50 and no threshold contributions. While there is no 
effect on the value of as{mz) the presence of the states in the 50 and 50 representations affects 
the value of the unified coupling at Mqut and also spoils the asymptotic freedom of SU(5) 
beyond Mqut- We find it suggestive that the solution of the doublet triplet splitting problem 
in terms of the 50, 50 and 75 representations automatically leads to improve the prediction of 
Oisi'^^z) and at the same time relaxes the constraints from proton decay. We now discuss this 
issue in more detail. 

Defining in the MS scheme a = aQEoijnz)) sin^ 6* = siv? 0{mz), = as{mz), = 
a^lMcur) , from the one-loop renormalization group evolution of the couplings in the MSSM 



6 



with only two light Higgs doublets we have: 



^ 60 sin^ e-12 

^(0) ^ 28a 

^ 36 sin^ ^ - 3 

/ , ^(Q) , N 7r(3 — 8 sin^ 6) , , 

log {M^SuT/mz) = (16) 

where by a'-"^ we mean the quantity a in the leading log approximation (LO). The LO results 
are the same as in minimal SUSY SU{^). For = 127.934 and sin^^ = 0.231 one obtains 
4°^ = 0.118 and M^auT = 2-1 10^*^ GeV. 

To go beyond the LO one must include two-loop effects in the running of gauge couplings, 
threshold effects at the scale msusY, close to the electroweak scale, and threshold effects at 
the large scale Mqut- The unification scale Mqut is not univocally defined beyond LO and 
here we choose to identify it with the mass of the SU (5) superheavy gauge bosons. The results 
derived in ref. |0 can be written in the following form: 

af 

as = 



(1 + af 5) 
r , 1 1 ^susY 3 niT 

5 = k + —\og — log— ^ . 17 

27r niz Svr M^{jj^ 

In this expression for 6 the logarithmic term with msusv comes from particles with masses 
near msusv- Similarly, the logarithmic term with ttlt arises from particles with masses of 
order itit ~ Mqut- The definition of mx is the mass of Higgs colour triplets in the minimal 
model or the effective mass, defined in eq. (|^, that, in the realistic model, plays the same 
role for proton decay |l^. The term indicated with k = k{2) + k{SUSY) + ki^McuT) 
contains the contribution of two-loop diagrams to the running couplings, A;(2), the threshold 
contribution of states near msusv , k{SUSY), and the threshold contribution from states near 
Mgut, k^Mcur)- The threshold contributions would vanish if all states had the mass msusv 
or Mgut so that only mass splittings contribute to k{SUSY) and k^Mcur)- The values of k{2) 
and of k{SUSY) are essentially the same in the minimal and the realistic model: typical values 
are k{2) = -0.733 and k{SUSY) = -0.510. The value for k{SUSY) = -0.510 corresponds to 
the representative spectrum displayed in Table 1. The value of k^Mcux) is practically zero for 
the 24 of the minimal model while we have k{MGUT) = 1-857 for the 75 of the realistic model. 
Thus we obtain: 

k = -0.733 - 0.510 = -1.243 minimal model 
k = -0.733 - 0.510 + 1.857 = 0.614 realistic model (18) 

This difference is very important and makes the comparison with experiment of the predicted 
value of as{mz) much more favorable in the case of the realistic model. In fact for k large and 
negative as in the minimal model we need to take msusv as large as possible and m^ as small 
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as possible. But the smaller ttit the faster is proton decay. The best compromise is something 
like f^susY ~ 1 TeV and itlt ~ Mq^jj,, which leads to as{mz) ~ 0.13 which is still rather large 
and proton decay is dangerously fast. This is to be confronted with the case of the reahstic 
model where k is instead positive and large enough to drag CKs{mz) below the observed value. 
We now prefer my to be larger than Mqij'p by typically a factor of 20-30, which means a factor 
of 400-1000 of suppression for the proton decay rate with respect to the minimal model. For 
example, for msusY ~ 0.25 TeV and itlt ~ 6 ■ 10^^ GeV we obtain as{mz) ~ 0.116 which is 
acceptable. The predictions for as versus {msusY,^T) in the minimal and the missing doublet 
models are shown in fig. 1. 



sparticle 


mass^ 


gluinos 


(2.7mi/2)' 


winos 


(0.8mi/2)2 


higgsinos 


/^^ 


extra Higgses 




squarks 


ml + 6m?/2 


(sleptons) L 


ml + 0.5m^/2 


(sleptons) 


ml + 0.15mf/2 



Table 1. Representative SUSY spectrum. SU{2) ® U{1) breaking effects arc neglected. 
The additional freedom related to the parameters mo, mi/2, /i, mn is here fixed by choosing 
O.Smo = 0.8m 1/2 = 2/x = mn and taking as a definition msusY = ^h, so that all particle 
masses can be expressed in term of msusY- This parametrization leads to k{SUSY) = —0.510. 

Clearly, there is an uncertainty on k{SUSY), related to the possibility of varying both the 
parametrization used and the relations assumed between the parameters mo, mi/2, A'-, ^h- An 
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Figure 1: Contours of 03 = as{mz) in the plane (m5[/5y,mr), for the minimal SU{b) and the 
missing doublet model. The SUSY spectrum is parametrized as in Table 1. 
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overall uncertainty of ±0.005 on 03 has been estimated in ref. |T9|| by scanning many models 



consistent with electroweak symmetry breaking, a neutral lightest supersymmetric particle and 
sparticle masses above the experimental bound and below ~ 2 TeV ^. 

Another source of uncertainty on is related to the unknown physics above the cut-off 
scale of the grand unified theory. There can be threshold effects due to new heavy particles or 
even non perturbative effects that arise in the underlying fundamental theory. These effects can 
be estimated from the non-renormalizable operators, suppressed by (y)/A, that split the gauge 



couplings at the scale Mq^j, |2^. For generic, order one coefficients and barring cancellations 
among different terms, the presence of these operators may affect 03 by additional contributions 
of about ±0.005Mpj/A [^. As we will see the model under discussion requires A ^ O.lMp;. 
Therefore, to maintain the good agreement between the experimental and predicted values of 
03, we need a suppression of this contribution by about a factor of 10, which we do not consider 
too unnatural. 

In view of the large theoretical uncertainties on we cannot firmly conclude that the gauge 
coupling unification fails in the minimal model while it is completely successful in the realistic 
one. However we find very encouraging that, by solving the doublet-triplet splitting problem 
within the missing partner model, acceptable values of 03 can be easily obtained and that they 
are directly related to a proton lifetime potentially larger than in the minimal model. 

Due to the large matter content the model is not asymptotically free and the gauge coupling 
constant blows up at the scale 

Taking into account that threshold effects modify 05 with respect the LO value af^ ^ 1/24, 
we find that the pole occurs near 10^^ GeV, the precise value depending on the details of the 
heavy spectrum. 



4 Yukawa couplings 

To reproduce the fermion mass spectrum we must further specialize the Yukawa couplings 
introduced in eq. (||): 

^ We stress that in the present analysis we are assuming universal soft breaking parameters at the cut-off 
scale. If we relax this assumption, then a larger range for as could be obtained. Indeed 03 is particularly 
sensitive to wino and gluino masses, that, at the electroweak scale, are approximately in the ratio 1:3 when a 
universal boundary condition on gaugino masses is imposed. It has been observed that, by inverting this ratio, 
a negative contribution of about —0.01 to 03 is obtained |2l|. 
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+ ^iG, ^5a^"-^^^lG'M^l + ... (20) 

where (r = (i, t/, M, 50) is proportional to Gr{{X),Q) of eq. (|). We have exphcitly 
introduced a term hnear in y, whose couphngs are described by the (X)-dependent matrix Fd. 
There are other terms hnear in F/A, not exphcitly given above. In particular we may insert 
y/A in the renormalizable term providing masses to the up type quarks. We neglect such a 
term since, on the one hand, the matrix is already sufficient to correctly reproduce the up 
quark masses and, on the other hand, this operator would not significantly modify the results 
for proton decay. As we will show, Y/K is close to 0.1 in our model and higher order terms 
in the Y/ A. expansion can be safely neglected. The interaction term involving H-^^f^^^ does not 
contribute to the fermion spectrum, but it will be relevant for the proton decay amplitudes. 

The term linear in Y in the previous equation is sufficient to differentiate the spectra in the 
charged lepton and down quark sectors. We get the following Dirac mass matrices: 

mu,u = yu,u , rnd,e = yd,e ' (2^) 



Vu = Gu yu = G^ 

(Y) (Y) (22) 

yd = Gd+^-T^Fd y^ = Gd-3^—^Fd , 



where f„, Vd and (Y) parametrize the vevs of H, H and Y respectively. The fermion spectrum 
can be easily fitted by appropriately choosing the numerical values of the matrices G„, G^, Fd, 
Gy and Gm- The most general fitting procedure would leave a large number of free parameters. 
Here we limit ourselves to the discussion of one particular example. In agreement with eq. (|12]), 
we take: 
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G 



M 



(1.50 + 0.55i)A2 (1.41 + 1.19i) (0.35 - 1.53i)A 
(1.41 + 1.19i) 

(0.35 - 1.53i)A 1.26 + 1.48i 



(27) 



While the generic pattern of G„, Gd, Fa, G^ and Gm is dictated by the U{1) flavour 
symmetry, the precise values of the coefficients multiplying the powers of A are chosen to 
reproduce the data. We take A = (X)/A = 0.25, tan/? = 1.5 and M = 0.9 ■ 10^^ GeV. In SU(5) 



the matrix G„ contains two additional phases [23|, 0i and 02 that have been set to zero in eq. 



These phases do not affect the fermion spectrum but enter the proton decay amplitude. 
When discussing the proton decay we will analyze also the dependence on 0i and 02- 

From the above matrices we obtain, at the unification scale: 

mt = 200 GeV = 0.27 GeV m„ = 0.9 MeV 

nib = 1-0 GeV m, = 26 MeV ma = 1.1 MeV (28) 
= 1.1 GeV = 71 MeV = 0.34 MeV 

|K«|=0.22 |K6|= 0.0022 11^61= 0.052 J = 1.9 • 10^^ , (29) 

where J is the CP- violating Jarlskog invariant. In the neutrino sector, we find: 

mi = 0.81 ■ IQ-hV , m2 = 0.88 • lO'W , mg = 0.061 eV . (30) 

More precisely: 

Ami, = ml- ml = 1.1 ■ 10~' , Aml,^ = ml - ml = 3.7 ■ 10'^ . (31) 
The neutrino mixing angles are 

^12-^ , ^23-^ , ^13 = 0.06 . (32) 

Until now we have not specified the value of {Y)/A. We know that {Y) should be around Mq^j.. 
The cut-off A cannot be too close to (Y), otherwise most of the spectrum of the model would 
lie beyond the cut-off. At the same time A cannot be too large: it is bounded from above by 
the scale at which the SU{b) gauge coupling blows up, which, as we see from eq. (|^), occurs 
more or less one order of magnitude below the Planck mass. This is welcome. If, for instance, 
we take {Y)/A = 0.05 — 0.1, it is reasonable to neglect multiple insertions of Y in the Yukawa 
operators. At the same time, in the example given in eq. (|25|) , the coefficients of the powers of A 
in Fd remain of order one, even if {Y)/A is as small as 0.05. A too large cut-off would have been 
ineffective in separating down quarks from charged leptons unless we had chosen unnaturally 
large coefficients in the allowed operators. What is usually considered a bad feature of the 
missing partner mechanism - the lack of perturbativity before the Planck mass - turns out here 
to be an advantage to provide a correct description of the fermion spectrum. 



The neutrino sector is quite similar to one of the two options described in ref. |T^. We 
obtain a bimaximal neutrino mixing with the so-called LOW solution to the solar neutrino 
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problem 0. Within the same U{1) flavour symmetry considered here we could as well reproduce 
the vacuum oscillation solution, by appropriately tuning the order one coefficients in G^, and 
Gm- We recall that the value of M required to fit the observed atmospheric oscillations is 
probably somewhat small in the context of SU{5), where a larger scale, closer to the cut-off A, 
is expected. This feature might be improved by embedding the model in SO (10) where M is 
directly related to the B — L breaking scale. 

In conclusion, the known fermion spectrum can be reproduced starting from a superpotential 
with order one dimensionless coefficients. Mass matrix elements for charged leptons and down 
quarks match only within 10 — 20%, due to {Y)/A ?s 0.1, and this produces the required 
difference between the two sectors. The neutrino mixing is necessarily bimaximal in our model, 
with either the LOW or the vacuum oscillation solution to the solar neutrino problem. 



5 Proton Decay 



Similarly to the case of minimal SU{5), we expect that the main contribution to proton decay 
comes from the dimension five operators ||2^ originating at the grand unified scale when the 

. We denote the colour triplets contained 



colour triplet superfields are integrated out p5 

H3d 



in H, H, H,o, by H- 



H'^^ and H'^^, respectively. The part of the superpotential 



depending on these superfields reads: 



w 



+ 

+ 
+ 



—-QGuQ 



U'G^.E' 



m 



3u 



QCL + WDD' 



--QG-,Q + WG-,E^ 



where 



C 



-Ga 



D = G„ 



(33) 



(34) 



A ' A 
and Q, L, W^, and E'^ denote as usual the chiral multiplets associated to the three fermion 
generations. Notice that, at variance with the minimal SU{'o) model, an additional interaction 
term depending on H'^^^ is present. By integrating out the colour triplets we obtain the following 
effective superpotential: 

1 

irtrp 



w, 



eff 



QAQ QGL + WBE" U^DD' 
where hit has been defined in eq. (^, 

Gu 



+ 



(35) 



B 



-2A 



(36) 



^Latest preliminary results from Super-Kamiokande including constraints from day-night spectra, seem 
in fact to prefer bimaximal neutrino mixing and to revamp interest in the LOW solution. 
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and dots stand for terms that do not violate baryon or lepton number. Minimal SU{5) is 
recovered by setting G-^ = = 0. In that case the matrices A, B, C and D are determined 
by Gu (in which now also (pi and 02 play a role) and Gd and therefore strictly related to the 
fermionic spectrum . In our case we have a distortion due to the terms proportional to G^ 
and Fd- These distortions have different physical origins. On the one hand the terms containing 
Fd are required to avoid the rigid mass relation of minimal SU{5). They are suppressed by 

On the other hand the terms containing 



{Y)/A and we expect a mild effect from them p8 



G50 are a consequence of the missing doublet mechanism and they might be as important as 
those of the minimal model. 



At lower scales the dimension five operators give rise to the four fermion operators relevant to 



proton decay, via a "dressing" mainly due to chargino exchange p4|, |2^, ^ . When considering 
the operators QQQL the main contribution, here called Li, comes from the exchange of a wino 
Q. Charged higgsino exchange provides instead the most important dressing of the operators 
V^V^D^E^ We term L2 the leading four-fermion operator in this case. Beyond Li and L2 
other 6 operators are generated by chargino exchange, 3 from QQQL and 3 from U'^U'^D'^E'^. 
The contributions to the proton decay amplitudes from these operators are suppressed by at 
least with respect to those associated to Li and L2, and can be safely neglected in the 
present estimate. Li and L2 are given by: 



L2 



Cljki i^idi) (djUk) + h.c. 

^ijkl {diUj] 



+ h.c. 



where i, j, /c, I are generation indices and 



(37) 



^ijkl 



^2 



iriT Svr 



[/i,(M, d) + /^,(m, e)] Pij Qki 

2 fhi'^'^y '^kl 



(3^ 



P 



2 Lji Ld 



Q 



Lie Le 



■'u,d,e 



S = L\ylBylLl T = Rl 1) R\ . 
and Ru,d,e are the unitary matrices that diagonalize the fermion mass matrices: 



(39) 



LlvuRi = iVu] 



diag 



'^dVdR, 



(y. 



djdiag 



M = (y, 



ejdiag 



(40) 



with {yf)diag (/ = u,d,e) diagonal and positive. The quark mixing matrix is Vckm = L\^Ld. 
K12 are constants accounting for the renormalization of the operators from the grand unification 
scale down to 1 GeV |2^, In our estimates we take Ki = K2 = 10. Finally, / is a 

function coming from the loop integration: 



/c(a, h) 



rrir 



2 

mi 



log 



2 

mi 

b 



m% 



m% 



log 



2 

mi 

b_ 

2 

mi 



(41) 



a c c b c c _ 

^Gluino dressing contributions cancel among each other in case of degeneracy between first two generations 
of squarks Q 
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We parametrize the partial rates according to the resuhs of a chiral lagrangian computation 
[p^]. The rates for the dominant channels are given in Table 2. 



To estimate the proton decay rates we should specify some important parameters. First of 
all the mass rriT- From the discussion of the threshold corrections we know that a large value 
for rriT is preferred in our model. We should however check that this value can be obtained 
with reasonable choices of the parameters at our disposal. The unification conditions constrain 
(Y) in a small range around 10^^ GeV 0. Then the cut-off scale and the X vev are fixed by the 
phenomenological requirements (Y) ^ 0.05 A and A = {X)/A = 0.25. Large values of rriT could 
be obtained either by taking 711x^711x2 large or by choosing a small m^. This last possibility is 
however not practicable, since the gauge coupling blows up at approximately 20 m^: it is 
not reasonable to push below 5 ■ 10^^ GeV. 



channel 


rate 




Xk 


/3 


&Ch,, + {l + i^{D + 3F))c^,,, 




2 


K+Vr 


Xk 




/3 










2 
3 


^DC!,,, + {l + i^{D + 3F))ci,,, 


2 






X^\p{l + D + F)Cl^^^\^ 




Xtt 


p{l + D + F) Cl,,^ + a{l + D + F) C^,,,, 


2 



Table 2. Proton decay rates. We define Xx^n = ijnl — '"^|',7r)^/(327rmp/^); mp, rriK, ttIt, are 
the proton, and 7r+ masses; is an average baryon mass, /tt is the pion decay constant; 
D and F are coupling constants between baryons and mesons in the relevant chiral lagrangian; 
(3 and a parametrize the hadronic matrix element. In our estimates we take mp = 0.938 GeV, 
7nB = 1.150 GeV, m^^ = 0.494 GeV, = 0.140 GeV, = 0.139 GeV, D = 0.8, F = 0.45 
3|], /5 = -a = 0.003 GeV=^ |6|. 



The only remaining freedom to obtain the desired large value for 7nx is represented by the 
coefficients C2 and C3 that, however, cannot be taken arbitrarily large. We should also check 



^More precisely, after the inclusion of threshold corrections and two loop effects, the last equality in eq. (16) 



becomes a relation among the masses of the super-heavy particles and the leading order quantity Mq^j.. This 
relation limits the allowed range for the heavy gauge vector bosons and indirectly pushes the 5[/(5) breaking 
vev (F) close to lO^^. 
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my <Y> <X> A 111^ ^pv,nck 

^5 ^ 4^ ^ ^ ^ — ^ — 

10 10' Mqut 10" 10'^ (GeV) 

Figure 2: Heavy sector of the spectrum. 

that all the heavy spectrum remains below A and this requirement imposes a further constraint 
on our parameters. A choice that respects all these requirements is provided by []: 

ci = 0.004 , C2 = C3 = 19.3, C4 = 0.7 , (42) 

(X) = 3.0 ■ 10^^ GeV, (Y) = 5.7 ■ 10^^ GeV (A = 1.2 ■ 10^^ GeV) , (43) 
which leads to: 

Mgut = 2.9 ■ 10^*^ GeV , = 2.0 ■ 10^*^ GeV , (44) 

rriT, = 1.2 • 10^^ GeV , rriT^ = 1.0 ■ 10^^ GeV , my = 6 ■ 10^^ GeV . (45) 

The heavy sector of the particle spectrum is displayed in fig. 2. With the above values we also 
obtain: 

^j^-5.5 . (46) 
C4 (X) 

To evaluate the loop function / we also need the spectrum of the supersymmetric particles. 
As an example we take here the same spectrum considered in table 1, with msusY = 250 GeV. 
This leads to a squark mass of about 800 GeV, a slepton mass of approximately 350 GeV a 
wino mass of 250 GeV and a charged higgsino mass of 125 GeV. 

The matrix G-^ is not directly related to any accessible observable quantity and the only 
constraint we have on it comes from the U{1) flavour symmetry that requires the following 
general pattern: 

■A^ A^" 

A6 A^ A3 . (47) 
A^ A3 A 



^50 



The texture for G-^ has an overall suppression factor A compared to Therefore the 
contribution of G^ to the matrices A and B is comparable or even slightly larger than the 
minimal contribution provided by Gu- The interference between the amplitude with Gu and 
the one with G-^ can be either constructive or destructive, depending on the relative phases 
between the two terms. We have scanned several examples for G-^, obtained by generating 
random coefficients for the order one variables in eq. (0). By keeping fixed all the remaining 
parameters we obtain a proton decay rate in the range 8 ■ 10^^ — 3 ■ 10^^ ys for the channel 



'^It is clear that the extreme values of the coefficients C2 and C3 here adopted raise doubts on the validity of 
the perturbative approach that has been exploited in several aspects of the present analysis. We adhere to this 
extreme choice also to show the difficulty met to obtain acceptable phcnomcnological results within a not too 
complicated scheme. 
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p — >■ K~^V and a rate between 2-10^^ ys and 8-10^^ ys for the channel p — > vr+F. For comparison, 
considering the same choice of parameters but setting G-^ = Fd = 0, we obtain 9 ■ 10^^ ys and 
2- 10'^^ ys respectively, for the above channels. The present 90% CL bound on t/ BR[p —>■ K^v) 
is 1.9 10^^ ys 0. These estimates have been obtained by setting to zero the two physical phases 
01, 02 contained in the matrix These additional parameters may increase the uncertainty 
on the proton lifetime. For instance, in minimal SU(5), the proton decay rates for the channels 
considered above, change by about one order of magnitude when 0i and 02 are freely varied 
between and 27r. Even when the inverse decay rates for the channels K^v and tt^v are as 
large as 10^^ ys, they remain the dominant contribution to the proton lifetime. Indeed, since 
the heavy vector boson mass, Mqut, is equal to 2.9 10^^ GeV in our model, the dimension 6 
operators provide an inverse decay rate for the channel e^vr^ larger than 10^^ ys. 

The effective theory considered here breaks down at the cut-off scale A. We expect additional 
non-renormalizable operators contributing to proton decay amplitudes from the physics above 
the cut-off. By assuming dimensionless coupling constants of order one, in unified models 
without flavour symmetries the proton lifetime induced by these operators is unacceptably 
short, even when A = Mpi ^ |7|- In our case these contributions are adequately suppressed 



by the U{1) symmetry. If we compare the amplitude Anr induced by the new non-renormalizable 
operators with the amplitude ^3 coming from the triplet exchange, we obtain, for the generic 
decay channel, 

^«£HSLf<I>y«i . (48) 

A3 \ A J 

This supports the conclusion that a proton lifetime range considerably larger than the one 
estimated in minimal models is expected in our case. 



6 Conclusions 

We have constructed an example of SUSY SU{5) GUT model, with an additional U{1) flavour 
symmetry, which is not plagued by the need of large amounts of fine tunings, like those as- 
sociated with doublet-triplet splitting in the minimal model, and leads to an acceptable phe- 
nomenology. This includes coupling unification with a value of as{mz) in much better agree- 
ment with the data than in the minimal version, an acceptable pattern for fermion masses 
and mixing angles, also including neutrino masses and mixings, and the possibility of a slower 
proton decay than in the minimal version, compatible with the present limits (in particular the 
limit from Super-Kamiokande of about 2 • 10^^ ys for the channel p K^u). In the neutrino 
sector the present model is a special case of the class of theories discussed in ref. |T^. The 



preferred solution in this case is one with nearly maximal mixing both for atmospheric and 
solar neutrinos. The U{1) flavour symmetry plays a crucial role by protecting the light doublet 
Higgs mass from receiving large mass contributions from higher dimension operators and by 
determining the observed hierarchy of fermion masses and mixings. Of course, the U{1) symme- 
try can only reproduce the order of magnitude of masses and mixings, while more quantitative 
relations among masses and mixings can only arise from a non abelian flavour symmetry. 
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A remarkable feature of the model is that the presence of the representations 50, 50 and 75, 
demanded by the missing partner mechanism for the solution of the doublet-triplet splitting 
problem, directly produces, through threshold corrections at Mqut from the 75, a decrease of 
the value of as{mz) that corresponds to coupling unification and an increase of the effective 
mass that mediates proton decay by a factor of typically 20-30. As a consequence the value 
of the strong coupling is in better agreement with the experimental value and the proton 
decay rate is smaller by a factor 400-1000 than in the minimal model. The presence of these 
large representations also has the consequence that the asymptotic freedom of SU{5) is spoiled 
and the associated gauge coupling becomes non perturbative below Mpi. We argue that this 
property far from being unacceptable can actually be useful to obtain better results for fermion 
masses and proton decay. 

Clearly such a model is not unique: our version is the simplest realistic model that we 
could construct. We think it is interesting because it proves that a SUSY SU{5) GUT is not 
excluded and offers a benchmark for comparison with experiment. For example, even including 
all possible uncertainties, it is difficult in this class of models to avoid the conclusion that 
proton decay must occur with a rate which is only a factor 10-50 from the present bounds. 
Failure to observe such a signal would require some additional specific mechanism in order to 
further suppress the decay rate. Finally it is a generic feature of realistic models that the region 
between Mqut and Mpi becomes populated by many states with different thresholds and also 
non perturbative phenomena occur. This suggests that the reality can be more complicated 
than the neat separation between the GUT and the string regime which is postulated in the 
simplest toy models of GUT's. 
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